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QCD 相 变 的 戴 森 - 施 温 格 方程 方法 研究 
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摘要 ”量子 色 动 力学 (Quantum Chromodynamics,QCD) 在 有 限 温度 有 限 密度 平面 内 具有 丰富 的 相 结构 ,关于 相 
边界 曲线 .临界 终点 (CCritical Endpoint,CEP) 的 位 置 以 及 各 相 的 热力 学 性 质 、 状 态 方程 等 的 研究 是 目前 理论 和 实 
仿 上 的 重点 问题 ,寻找 QCD 相 变 信号 ,特别 是 确定 CEP 的 位 置 ,成 为 了 相对 论 性 重 离子 对 撞 机 当前 和 未 来 实验 
计划 的 主要 目标 。 本 文 简要 小 结 在 QCD 的 戴 森 - 施 温 格 方程 (Dyson-Schwinger equation, DSE) 方 法 框架 下 , 关 
于 手 征 对 称 性 恢复 相 变 和 退 禁 闭 相 变 研究 的 进展 和 主要 结果 ,着 重 介绍 关于 DSE 的 截断 方案 的 探索 和 完善 、 
QCD 相 变 的 手 征 磁化 率 判 据 和 推广 的 施 温 格 (Schwinger) 函 数 判 据 、 及 由 之 确定 的 手 征 和 禁闭 两 种 相 变 的 相 
图 。DSE 方 法 预言 的 手 征 相 变 的 CEP 位 置 大 致 位 于 (7,w,) (110,600) MeV ,与 泛 函 重 整 化 群 等 方法 预言 的 结 
果 一 致 。 这 些 结果 对 于 重 离子 对 撞 实 验 搜寻 CEP 信 号 具有 指导 价值 。 

关键 词 QCD 相 变 ， 手 征 对 称 性 动力 学 破 缺 ， 解 禁 财 ， 戴 森 - 施 温 格 方程 ， 手 征 磁化 率 
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QCD phase transitions using the QCD Dyson-Schwinger equation approach 
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2(School of Physics, Peking University, Beijing 100871, China) 


Abstract The use of the relativistic heavy ion collision experiment has extended our insights into the diverse 
possibilities available to a truly strongly-interacting system. The main goal of this experiment is to describe the 
properties of the different phases of quantum chromodynamics (QCD) and to chart the QCD phase diagram on the T- 
mu plane. For the phase diagram, apart from the general phase boundary lines, some specific characteristics such as 
the possible critical endpoint (CEP), associated coexistence region, and strongly-coupling quark-gluon plasma 
(sQGP) have to be identified. Here, the CEP separates the first-order phase transition from the second-order transition 
(or crossover) when the case beyond the chiral limit is considered. However, convincing signals have not yet been 
obtained using the relativistic heavy ion collider (RHIC) experiment. Theoretically, strong interaction systems hold 
significant features: asymptotic freedom in the ultraviolet region, dynamical chiral symmetry breaking, and 
confinement in the infrared region. Such features can be uniformly displayed in the phase structure of the matter in 
the temperature 7 and chemical potential planes. Consequently, several investigations have been experimentally and 
theoretically performed. However, the strong coupling feature in the low-energy region prevents the use of 
perturbative calculation methods, which creates the need for the development of nonperturbative approaches. 


Additionally, lattice QCD simulations have been widely implemented; however, the "sign problem" delays the 
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progress in the large chemical potential region. Therefore, the Dyson-Schwinger equation (DSE) equation method 
and functional renormalization group approach, which inherently include both dynamical chiral symmetry breaking 
(DCSB) and confinement, play an important role. The QCD DSE approach is a method based on the continuum 
quantum field theory. The new criteria were proposed based on the DSE and studied using the deconfinement and 
Chiral symmetry restoration phase transition of QCD. Currently, functional methods can be used to provide a reliable 
estimation of the CEP location. First, reliability is achieved using a thorough investigation of the truncation of the 
DSE, state of the art truncation is then performed causing a converging result between the different methods, and the 
predication of the lattice simulation at low chemical potential is confirmed. The results show a fast convergence of 
the truncation owing to the infrared fixed point of the QCD coupling, which allows the capturing of the QCD running 
behavior using a finite set of two- and three-point Green functions. The estimated location of the CEP based on the 
current computation is 4p at 600~650 MeV and T at 100~110 MeV. The existing functional QCD methods are non- 
perturbative continuum methods that are capable of simultaneously describing both the DCSB and confinement. 
Although they are limited by the truncations, the use of functional QCD approaches has resulted in progress in the 
study of the QCD phase structure and thermal properties, where a complete phase diagram and related thermal 
properties have been obtained in a large chemical potential range, which can provide a reference for the exploration 
of the QCD features. Most of the theoretical studies using effective models or certain truncations have observed the 
existence of the CEP; however, the determination of its location is still a work in progress because it varies based on 
the computation. Moreover, searching for QCD phase transition signals, particularly the CEP, is the main goals of 
current and future experimental programs on the relativistic heavy ion collider. 

Key words QCD phase transition, Dynamical chiral symmetry breaking, Deconfinement, Dyson-Schwinger 


equations, Chiral susceptibility 


物理 学 “ 见 物 讲理 、 依 理 造物 ”的 学 科 性 质 和  ” 相 变 及 有 限 温 度 密度 下 各 相 物 质 的 状态 和 性 质 将 深 


范畴 知道 ,对 我 们 所 处 宇 害 中 的 可 见 物质 及 其 质量 ”化 我 们 关于 强 相互 作用 基本 规律 的 认识 ,并 揭示 可 
的 产生 和 演化 的 研究 自然 是 物理 学 中 最 重要 的 课 见 物质 及 其 质量 起 源 的 奥秘 ””。 


题 。 


按照 现在 的 认识 ,我 们 现在 所 处 的 宇宙 起 源 于 由 于 QCD 相 变 是 发 生 在 QCD 的 非 微 扰 能 区 的 


130 多 亿 年 之 前 的 一 次 大 爆炸 ,经 过 平稳 演化 阶段 ”夸克 胶 子 物质 的 状态 演化 现象 ,因此 ,人 们 利用 相对 
之 后 出 现 暴 涨 。 相 应 于 暴涨 引起 的 温度 降低 和 空间 论 性 重 离子 对 撞 对 之 进行 实验 研究 、 利 用 非 微 扰 
脱 胀 ,不 可 分 辨 的 夸克 和 电子 等 轻 子 出 现 ;再 经 降 QCD 方法 对 之 进行 理论 研究 。 目 前 ,理论 上 发 展 起 
温 , 夸 殉 因 著 得 不 同 的 流质 量 而 可 以 区 分 ;再 进一步 。 来 的 非 微 扰 QCD 方法 有 离散 场 论 层次 上 的 格 点 


降温 ,可 见 物质 的 主要 组 成 成 分 


强 子 出 现 ,相应 ” QCD 模 拟 方法 、 连 续 场 论 层 次 上 的 戴 森 - 施 温 格 方程 


的 可 见 物质 具有 了 质量 。 显 然 , 可 见 物 质 的 最 底层 (Dyson-Schwinger Equation,DSE) 方 法 和 泛 函 重 整 
的 组 分 粒子 是 夸克 和 胶 子 ,其 产生 与 演化 过 程 即 夸 化 群 (Functional Renormalization Group, 人 RG) 方 法 
克 ( 胶 子 ) 物 质 的 状态 的 演化 过 程 。 与 统计 物理 学 原 ” 及 多 种 唯 象 模型 方法 。 经 过 50 年 的 努力 ,理论 和 实 
理 对 应 ,这 一 过 程 即 强 相 互 作用 物质 ( 亦 常 简称 核 物 验 对 于 有 限 温 度 密 度 的 QCD 相 变 和 相 图 的 研究 都 
质 ) 的 相 变 。 由 于 表述 强 相互 作用 的 基本 规律 的 理 ”取得 了 很 大 的 进展 上 9。 一 般 而 言 ,在 温度 密度 平 
论 是 量子 色 动 力学 (Quantum Chromodynamics， 面 上 ,QCD 相 变 分 为 两 种 : 手 征 对 称 性 破 缺 与 手 征 


QCD) ,因此 ,该 过 程 常 被 称 为 QCD 相 变 。 由 目前 认 对 称 性 恢复 之 间 的 相 变 ( 常 简称 为 手 征 相 变 )、 禁 闭 
识 的 强 子 的 线 度 和 强 子 物理 的 基本 原理 知 ,发 生 ”与 退 禁 闭 之 间 的 相 变 (禁闭 相 变 )。 两 种 相 变 在 不 同 


YE 


QCD 相 变 的 能 标 在 10° MeV ,也 就 是 在 QCD 的 典型 ” 的 判 据 下 表现 并 不 相同 。 在 零 化 学 势 下 的 格 点 
的 非 微 扰 能 区 ,因此 ,抽象 来 讲 , 早 期 宇宙 强 相 互 作 = QCD 模拟 计算 ”中 ,使 用 Polyakov 圈 判 据 发 现 两 
用 物质 的 演化 乃 强 相互 作用 物质 在 能 标 跑 动 下 展现 。 种 相 变 具有 较 好 的 一 致 性 ,而 在 DSE 的 计算 却 给 出 
出 截然 不 同 的 自由 度 涌现 出 丰富 而 复杂 的 现象 的 Polyakov 圈 在 有 限 密度 区 域 退 禁闭 相 变 较 手 征 相 变 
表现 。 较 具体 地 ,这 一 演化 即 强 相 互 作用 物质 随 着 。 温度 仍 会 有 一 定 差别 的 结论 , 退 禁 闭 温度 略 高 于 
温度 密度 变化 由 不 禁闭 不 可 见 的 具有 和 手 征 对 称 性 的 。 手 征 相 变温 度 ,出 现 “quarkyonic” 相 ,这 与 大 N. 极 限 


夸 殉 胶 子 相 到 色 禁 闭 的 手 征 对 称 


见 的 强 子 相 的 相 变 过 程 ”。 由 此 可 知 ,研究 QCD 。 子 谱 函数 正定 性 的 分 析 则 得 到 轻 味 系 统 的 手 和 有 
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高 飞 等 : QCD 相 变 的 戴 森 - 施 温 格 方程 方法 研究 


与 退 禁 闭 相 变 一 致 的 结论 "下 ,关于 关联 函数 的 计 
算 和 分 析 也 说 明 退 禁闭 温度 相 变 温度 与 手 征 相 变温 
度 相同 与 否 依赖 于 系统 的 味道 ”。 

考虑 到 对 手 征 相 变 的 定义 及 性 质 的 描述 非常 清 
晰 并 且 是 QCD 相 变 的 主要 模式 ,因此 人 们 对 手 征 相 
变 的 相 图 已 进行 了 大 量 的 理论 研究 下 。 在 高 温 低 
化 学 势 情况 下 , 格 点 QCD 模 拟 和 DSE 方 法 等 的 计算 
都 表明 以 温度 为 主要 因素 驱动 的 QCD 的 手 征 相 变 
为 连续 过 渡 ”" 。 然 而 ,由 于 符号 问题 , 格 点 QCD 
很 难 探索 大 化 学 势 的 区 域 ,因此 关于 大 化 学 势 情 况 
下 的 临界 终点 (Critical Endpoint, CEP) 和 一 阶 相 变 
区 域 的 存在 性 的 研究 主要 是 通过 连续 场 论 层次 上 的 
DSEs "Al FRG" RZ A QCD 方法 来 实现 。 泛 函 
QCD 方法 是 不 仪 是 能 够 同时 描述 手 征 对 称 性 动力 
学 Wk tk (Dynamical Chiral Symmetry Breaking, 


1 1 2 7. 一 
£= (Fs) - 32 (0°43) + (iD - m)y + - 


a Ds Jes (1) 
HH: D, = 0, — igAit’, De = 0,6 - ig? f, Fa = 
0,43 — 0,42 + gf” ALAS. 

由 此 可 以 得 到 生成 泛 函 为 : 

Z ae fiue" + (ay + gn +J’ A +E + Ee) dy 

du = DyDyDA* D De (2) 
基于 对 生成 泛 函 的 泛 函 积分 内 对 某 场 的 变 分 为 
零 , 例 如: 


Ò iS + il (ny + gn +J:A)d'y _ 
[ax T° 0 G) 
即 得 相应 场 的 DS 方程 。 


DS 方 程 是 一 组 包含 无 穷 多 个 积分 方程 的 不 封 


DCSB) 和 禁闭 的 非 微 扰 连 续 方法 ,还 可 以 覆盖 大 范 
围 的 化 学 势 区 域 ,因此 在 研究 QCD 相 变 研究 方面 已 


闭 的 方程 组 ,它们 由 不 同 阶 的 格林 函数 耦合 而 成 。 
在 实际 应 用 中 ,人 们 需要 引入 一 定 的 截断 方案 才能 


取得 了 丰硕 的 进展 。 结 果 表 明 : 在 较 大 化 学 势 
下 ,以 化 学 势 为 主要 因素 驱动 的 相 变 是 一 级 相 变 ,在 
确定 相 边界 平衡 曲线 同时 ,这 些 计算 也 确认 了 连续 
过 渡 与 一 级 相 变 的 分 界 状态 , 即 CEP。 本 文 对 DSE 
方法 下 关于 QCD 相 变 研究 的 方案 、 核 心 技术 与 难 
点 ,以 及 取得 进展 等 予以 简要 介绍 。 


1 Dyson-Schwinger (DS) È fz SK QCD 
相 变 中 的 应 用 


1.1 DS 方程 的 一 般 形式 


Dyson-Schwinger(DS) 方 程 是 场 论 中 关于 格林 
函数 的 运动 方程 ,最 早 由 Dyson 和 Schwinger 建立 。 


对 DS 方程 进行 求解 。 原 则 上 ,所 有 的 观测 量 均 可 从 
格林 函数 中 提取 出 来 ,将 QCD 的 作用 量 代 入 ,并 求 
解 相 应 的 DS 方程 即 可 得 到 QCD 的 所 有 物理 性 质 。 
考虑 到 QCD 作用 量 的 具体 形式 ,n 点 的 格林 函数 的 
DS 方程 会 包含 (n+1) 点 和 (n+2) 点 的 格林 函数 ,而 且 
越 高 阶 的 格林 函数 的 DS 方程 形式 越 复杂 。 因 此 通 
常情 况 下 ,可 以 采取 一 些 合理 的 截断 方案 ,只 求解 最 
低 阶 格林 函数 的 DS 方程 。 为 了 得 到 关于 手 征 对 称 
性 破 缺 以 及 禁闭 的 信息 ,通常 用 到 的 是 最 低 阶 的 夸 
克 胶 子 传播 子 和 夸克 胶 子 顶 点 的 DS 方 程 \ 相 应 的 夸 
克 和 胶 子 以 及 顶点 的 费 曼 图 如 图 1 所 示 。 
为 了 研究 QCD Wa ZEA BRR din BE BE EP SAL 
还 需要 根据 有 限 温 度 场 论 的 原理 引入 温度 和 密度 自 
由 度 。 如 果 仅 考虑 平衡 态 的 信息 和 平衡 态 相 变 , 例 


通过 DS 方程 可 以 得 到 场 的 格林 函数 ,进而 得 到 所 有 
的 可 观测 量 。QCD 的 DS 方程 既 能 体现 QCD 的 手 
征 对 称 性 动力 学 破 缺 又 能 体现 色 禁 闭 ,是 一 种 具有 
很 好 的 普 适 性 和 QCD 基础 的 非 微 扰 连续 场 论 方法 ， 
在 对 强 相 互 作 用 系统 的 研究 中 取得 了 许多 有 意义 的 
结果 。 为 了 得 到 量子 场 的 DS 方程 ,首先 需要 给 出 场 
的 生成 泛 函 。 引 入 和 鬼 场 后 ,QCD 的 Faddeev-Popov 
拉 式 量 为 : 


b)~@r = + ew OrO@= 
y i 下 iD iD, iDoill iD 


WU A PS] DA Ae JE eh. BE E 2S JE n BS E EA OE 
IS, (Mg 22 S| A EI ve FEE A 
征 相 应 粒子 的 密度 。 对 于 系统 的 演化 信息 以 及 输 运 
性 质 等 非 平 衡 态 的 信息 , 则 需要 实时 温度 场 论 , 即 
Schwinger-Keldysh 形式 。 本 文 仅 介 绍 关 于 QCD 相 
图 的 研究 ,因此 局 限于 虚 时 温度 场 论 的 内 容 。 首 先 ， 
考虑 平衡 体系 的 算 符 4 的 平均 值 为 : 


eee 


图 1 SRT Se TAY DSE 77 TEN S AR as 


Fig.1 Feynman diagram for the DSE of the quark, gluon, and kuark-gluon interaction vertex 
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Tr (etA) 
其 中 :Z= Ty(e CAPRA TAE 
Eo onno 
Ta 
= f D[4] al f "agast Ao] (4) 
= [p[y] elit faxc a[o] 
这 表明 , 配 分 函数 相当 于 虚 时 空间 的 路 径 积 / 


虚 时 reL0,p]。 配 分 函数 中 的 求 迹 为 场 带 来 边界 条 
件 。 对 玻 色 子 , 有 格林 函数 : 


= S Tr(a) e Meng (5.0) ead 
(5) 


“Eee mr! B(x, t)6(¥p)| 
43.5 b(¥,B) = f(D, 0), Mt REIN AY 6 (x, r Ef 
IH} AEH AS Bh Ee Tal H S(p, o): 


周期 边界 条 件 得 :w, = 2nnT. 
对 于 费 米 子 , 最 后 一 步 交 换 歼 (x) 与 w(7) 会 产 
A HS Pw (9, B) = yy 0), FH: 


= ] dp sts inro 
vy (zr)= -5,| FIBO) 中 


=(2n + 1)aT (7) 
由 此 ,前 述 零 温 零 密 的 方法 依旧 适用 于 处 于 平 
衡 态 的 有 限 温 度 有 限 密度 体系 , 仅 O(4) 对 称 性 退化 
成 0(3) 对 称 性 ,相应 的 传播 子 以 及 顶点 的 形式 变 为 
0(3) 的 参数 化 形式 。 

1.2 ”夸克 胶 子 传 播 子 的 gap 方程 及 其 解 的 一 般 性 
质 与 相 变 判 据 
为 了 得 到 QCD 的 相 图 ,我 们 需要 对 夸克 胶 子 传 
播 子 以 及 夸克 胶 子 顶点 按照 洛 伦 芯 结 构 进 行 分 解 才 
能 数值 求解 。 其 中 ,夸克 传播 子 可 以 写 为 : 
S(p.e,) = 
lip - PA( p62) + B( B®) + ið C (Pa) 
(8) 


ag. 


相应 的 胶 子 传播 子 为 : 


SDI Po) = P*( Do.) Ar Po + P?( Po.) Acl Po, KO) 


其 中 :po: =(P, Q,), Q, = 2knT, PE PI I RÆ ô - 
双全 的 三 动量 和 能 量 分 量 : 
pP 
0, LV 三 4 
P? = Dp. 
mPa) = 5, PPL, yy = ij= 123 
P (10) 
PaP» 
P po) = bm = Ts = Pal Pa) 


为 具体 求解 夸克 传播 子 的 gap 方程, 我 们 需要 
对 夸克 胶 子 顶点 进行 一 定 的 截断 。 关 于 夸克 胶 子 顶 
点 的 截断 方案 决定 了 DSE 方 程 方 法 的 可 靠 性 ,也 是 
目前 DSE 方程 方法 研究 最 重要 的 部 分 。 从 洛 伦 效 
结构 来 说 ,顶点 横向 部 分 可 以 写 为 :T= 
(gp 如)TY, 其 中 Tr 是 一 组 完备 林 * 中 ,顶点 纵向 
部 分 有 4 项 ,可 由 瓦 得 -高 桥 恒等式 (Ward Takahashi 
Identity, WTD 完全 约束 为 Ball-Chiu 顶点 ,一 般 表 
IN: 

TE(G, onp, o,) = pros td ct 


1 1 
(p, + 1), 572. + gi), Ast xa pr + 1) Ac- in| 
其 中 : 


eae 1 
DAP 0P, o) = F(F 07 


Alg, onp, œo) = 


其 中 : FA, B, C, yl =y, yE YE =U, Yala U= 
(0,0,0,1) 

T:(tk) = s(t, - t- kk,) 
i 
A - t+ kk,) 
Ttk)= ky,- kk 
(Eu tin) tk - ki) 
iy , 
T= Ak) ake, 
T(tk)= t+ ky,- ket, 
Ti(t,k) = A k(t,- y,t) + gtd ki - tk) 
Ti(tk) = (kt, ~ ik,) + (ki = ik) 

ka = Qu Putu = Get Pa (11) 

目前 的 研究 表明 ,其 中 对 gap 方程 贡献 最 重要 
Me Tk) Ti, k) UK T(t, k) =, P 
T?(t, kE MRLE, TC, 大 ) 的 结构 来 自 于 
5) JI E A E AY DSc FS PARE CN) ,同时 也 被 横向 WTI 


<=> 


040015-4 


202306.00131v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


高 飞 等 : QCD 相 变 的 戴 森 - 施 温 格 方程 方法 研究 


的 研究 进一步 确认 后 ,对 于 gap 方程 的 贡献 也 尤其 
重要 。 目 前 来 讲 , 大 量 研究 基于 裸 顶 点 的 截断 方案 ， 
也 就 是 彩虹 -梯子 近似 方案 ,该 方案 可 以 保证 基本 的 
对 称 性 ,同时 便于 求解 。 然 而 ,由 于 缺乏 重要 的 
T(t, 大) 项 ,因此 简单 使 用 裸 顶 点 不 能 自然 地 得 到 手 
征 对 称 性 动力 学 破 缺 。 需 要 引入 相互 作用 模型 并 采 
用 较 大 的 耦合 常数 才能 在 真空 正确 的 层次 上 描述 厅 
标 和 矢量 介子 等 受 对 称 性 保护 的 强 子 的 行为 。 而 ! 
于 其 并 不 能 很 好 地 描述 相应 的 夸克 胶 子 场 的 跑 动 行 
为 ,因此 仅 能 定性 地 描述 QCD 相 图 。 目 前 基于 裸 顶 
点 的 相互 作用 模型 多 采用 Qin-Chang 模 型 : 
8? 8n°y, 


Gal’) = ph De” + 2 a- 
k 
lnljz 十 |1 + 
| a 


e Erni Yli? (12) 

其 中 :第 一 项 主要 表征 红外 相互 作用 ,耦合 强度 一 般 
1 

HN (Dw)? = 0.52 GeV, œ=0.5 GeV。 第 二 项 为 紫 

外 微 扰 项 ,其 中 t= e? - 1,m=0.5 GeV, y, = 12/25, 


Arey =0.234 GeV。 耦 合 强度 的 选取 由 真空 下 求解 


和 p 介 子 质量 以 及 衰变 常数 后 与 实验 值 对 比 来 
最 新 的 研究 表明 ,如果 同 时 考虑 
T(t, k) T(t, k) AR T(t, k) =A BRIE PR A 
夸克 胶 子 和 夸克 胶 子 顶点 DSE 时 ,QCD 的 跑 动 行为 
已 经 可 以 被 很 好 地 描述 。 因 此 ,目前 DSE 方 程 关于 
QCD 相 图 的 研究 已 经 可 以 定量 地 给 出 相 变 曲线 以 
K CEP 大 概 位 置 。 这 背后 的 深层 原因 是 QCD 耦合 
常数 的 红外 固定 点 使 得 相互 作用 顶点 处 于 近 共 形 极 
限 ,这 使 得 泛 函 QCD 方法 的 截断 方案 趋 于 稳定 , 保 
证 其 在 极 大 范围 能 标 下 均 可 以 准确 给 出 QCD 基本 
场 的 跑 动 行为 ,从 而 统一 自治 地 描述 强 相 互 作用 物 
AO, KMART UM QCD 的 有 效 荷 的 计算 中 得 
到 印证 ”。 图 2 给 出 了 关于 QCD 有 效 荷 的 计算 ,其 
中 实 线 和 折线 为 DSE 计 算 ”, 点 线 为 全 息 QCD 计算 
的 结果 2。 从 图 2 可 以 看 到 ,包含 两 点 和 三 点 格林 
函数 跑 动 行为 的 DSE 方 程 的 解 与 实验 观测 的 有 效 
荷 非 常 接 近 , 因 此 ,尽管 DSE 方 程 是 一 个 无 穷 耦 合 
的 方程 组 ,但 是 固定 点 的 标 度 现象 以 及 因此 带 来 的 
在 固定 点 附近 可 以 微 扰 展开 的 事实 使 得 高 阶 的 格林 
函数 对 最 终结 果 影 响 很 小 ,因此 ,DSE 方 法 的 截断 趋 
于 封闭 。 
在 截断 方案 确定 之 后 DSE 方 程 即 可 直接 数值 
求解 ,进而 给 出 手 征 对 称 性 动力 学 破 缺 以 及 色 禁 闭 


2Hall A / CLAS 
> JLab CLAS (2008) 


0 0.050.1 1 10 
K / GeV 


E2 ”过程 不 依赖 的 QCD A A” 


Fig.2 Process-independent effective charge"””! 


的 性 质 。 手 征 对 称 性 动力 学 破 缺 与 夸克 传播 子 直接 
相关 ,在 手 征 极限 下 , 当 夸 元 传播 子 的 标量 部 分 ,也 
就 是 质量 项 不 为 零 时 为 手 征 对 称 性 破 缺 的 Nambu 
相 , 质 量 项 为 零 时 为 手 征 对 称 的 Wigner 相 。 因 此 ， 
可 以 直接 地 定义 零 动量 的 质量 函数 8(p? = 0) 为 手 
征 极 限 下 的 序 参量 。 更 完整 的 序 参 量 定义 则 是 通过 
传播 子 求 迹 来 定义 手 征 凝 聚 : 


- (aa) = fusta) = mf eslo) 
(13) 
dtp B(p°) 


(2x) p?A’( p?) + Bp’) 
在 手 征 极限 下 , 手 征 凝 聚 是 恨 好 定义 的 。 然 而 ， 
在 超越 手 征 极限 时 , 式 中 的 积分 发 散 ， 


的 减 除 方案 。 首 先 , 利 用 求 和 规则 可 以 给 出 质量 
数 在 紫外 的 发散 行为 wm 
n BOPE EE 
= = 1 Ym 
Wilson on 7 GG rs 
L AA 
lin (Oaka) = 
因此 ,我 们 可 以 定义 超越 手 征 极限 的 序 参量 关 
Ai), 
(34) = (44), -mo (15) 


序 参量 可 以 通过 MEKAAR 势 给 出 相 变 
的 位 置 和 相 变 级 次 。 首 先 ,热力 学 势 可 以 展开 为 : 


T (aq) = TC) + aliay + EBay’ + 


zrla” (16) 
定义 磁化 率 为 ， 
~ aa) 


es, 


d= (17) 
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- (gq) Fe coe ~ Gay Fe eos. = Gay’ GB exc 


j at ra + 5y(aq)" 
_ 1 
` 
(49) es t G E ie a+ (Gay (SP) t- E 
PT 
7) or 
alga) man” 
因此 , 当 热 力学 势 二 阶 导 大 于 0, 也 就 是 系统 相 (po -r5 EP e op) 
处 于 热力 学 势 的 极 小 值 位 置 情况 下 ,对 应 于 手 征 磁 
化 率 大 于 0 的 相 为 稳定 相 ,对 应 于 磁化 率 小 于 0 的 相 j eae V 
为 不 稳定 相 。 两 相 共存 时 ,磁化 率 较 小 的 相 认 为 是 去 | do 0) enn 
亚 稳 相 。 据 此 ,利用 磁化 率 判 据 , 对 完全 非 微 扰 因而 ee 是 谱 函 数 正 定性 的 充 
无 法 确定 系统 热力 学 势 的 情况 ,我 们 也 可 以 给 出 其 。 分 条 件 ; 如 果 谱 函数 是 正定 的 , 则 Schwinger 函数 
手 征 相 变 的 相 边 界 曲 线 及 CEP 的 状态 。 A(z, p = 0) > 0; 如 果 出 现 禁闭 , 谱 函 数 存在 负 值 , 则 
关于 禁闭 ,目前 的 判 据 主 要 在 于 胶 子 传播 子 以 。 Schwinger 本 数 的 正定 性 会 被 破坏 。 这 一 判 据 可 以 


BRAK HY Polyakov # HIZ 吉 构 ,从 夸克 传播 子 判断 禁 进 步 f poten ， 直接 数值 求 导 可 以 看 到 Schwinger 


闭 的 主要 依据 是 人 克 谱 的 正定 性 ,首先 夸克 传播 了 。 函数 的 2n 阶 导数 的 正定 性 与 谱 函 数 的 正定 性 保持 


_ Fane p,(M?) - ipp,( M°) 1 rz 四 -(o+mr 
2s PFF one f dolo +u 00) 
o(p) = [PE LU) = ipM’) (21) 
o 2 p +M’? 因此 ， 我 们 可 以 通过 这 一 系列 导数 的 正定 性 来 

=0,( 三)+ 0,(p’) (19) 判断 谱 函 数 的 正定 性 ,从 而 得 到 夸克 禁闭 的 信息 


于 谱 函 数 的 信息 sy 需要 实时 温度 场 论 ,很 难 进 = > Zz FR ` FE 木 广 
行 非 微 扰 求 解 ,因此 ,人 们 通常 定义 Schwinger 函数 2? ”DS 方程 方法 下 的 关于 手 征 相 变 和 退 禁 闭 


来 判断 谱 函数 的 正定 性 鸣 。 在 有 限 温度 和 密度 下 ， ” 相 变 的 主要 结果 
可 以 通过 传播 子 的 谱 表 示 定 义 如 下 的 Schwinger 
函数 ; 21 手 征 相 变相 图 
首先 , 相 边 界 曲 线 可 以 从 热力 学 的 角度 进行 一 
般 的 判断 ,考虑 相 变 发 生 的 位 置 位 于 两 相 压 强 相同 
的 状态 , 即 有 P= Py WA: 
d(P, - P, far = (aPy/ðu = aP ,/ð A) au/aT + (aP,/aT -ðP „/ðT) = 0 (22) 
其 中 :压强 对 4 的 导数 即 为 粒子 数 密度 9Py94 = ny 行为 密切 相关 。 目 前 格 点 QCD 模拟 给 出 对 于 2+1 
xt TA FROME, MW aP T= sw。 一般 来 讲 , 手 。” 味 的 QCD 系统 ,x-0.015, 最 新 的 DSE 方 法 以 及 泛 函 
征 对 称 性 恢复 的 Wigner 相 的 精密 度 与 粒子 数 密度 BEARN ABA H T AWAR, k1 ai HA 


均 大 于 Nambu 相 的 相应 的 值 ,由 此 我 们 得 到 相 变 曲 “ ”不 同 QCD 方 法 给 出 的 最 新 结果 的 比较 。 


线 应 满足 如 下 关系 : 表 1， 相 变 曲线 曲率 
í 1 相 变 曲线 曲率 
W/oT = (Sy = Sy)/(My-Ny)<O (23) Tablet Curvature coefficients K 
这 表明 随 着 化 学 势 增 大 , 相 变 温度 会 降低 。 同 。 计算 方案 Methods Æ Coefficients SSS 
时 ,在 小 化 学 势 区 域 可 以 对 相 变 曲线 按 替 级 数 进行 DSE 0.015 0 
展开 : DSEE 9 0.023 8 
BoSS i : RG" 0.014 2(2) 
T( Ms)/T. 1 -K(M)/T, TAC uT.) (24) Lattice QCD"! 0.015 3(18) 
通过 上 述 分 析 可 以 看 到 ,x 应 为 正 值 。x 的 具体 Lattice QCD™ 0.014 4(26) 
数值 决定 相 变 曲 线 的 主要 行为 ,并 且 与 QCD 的 跑 动 Lattice QCD™! 0.015(451) 
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Chiral Susceptibility 


3 “一 些 不 同化 学 势 下 , 手 征 磁化 率 随 温度 变化 的 行为 中 
Fig.3 T dependence of the chiral susceptibility at different 
chemical potentials™ 

在 整个 温度 化 学 势 平面 内 数值 求解 夸克 传播 子 
的 gap 方程 后 可 以 直接 带 入 手 征 磁 化 率 的 相 变 判 所 
计算 其 随 温 度 和 化 学 势 ( 密 度 ) 变 化 的 行为 。 一 些 化 
学 势 下 磁化 率 随 温度 变化 的 行为 的 结果 如 图 3 中 展 


0.14 
从 这 长 = CEP ~ (0.11, 0.14) GeV| 一 7 smm 
-XWin 
0:10 Chiral Symmetry 
% 0.08 \ Restored 
9 0.06 Chiral Symmetry Broken \\ 
Sooo ` 
0.04 上 
0.02 Co-existence |; 


' (a) 
0 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 
uI GeV 


示 。 可 以 看 到 ,在 化 学 势 较 小 的 区 域 , 手 征 磁化 率 是 
一 个 光滑 函数 ,表明 QCD 相 变 实际 是 连续 过 渡 , 可 
以 认为 极 大 值 的 位 置 对 应 的 温度 即 硕 相 变 温度 。 随 
着 化 学 势 逐 渐 增 大 ,磁化 率 趋 于 发 散 ,发 散 的 位 置 对 
应 的 状态 即 为 CEP 的 状态 。 在 更 大 的 化 学 势 处 , 磁 
化 率 断 开 , 标 志 着 系统 进入 一 级 相 变 区 域 。 
首先 ,人 们 运用 QC 胶 子 模型 对 QCD 的 手 征 相 
变 进行 了 大 量 的 研究 。 由 于 截断 方案 的 限制 ,相关 
的 计算 仅 能 定性 确定 QCD 相 变 的 行为 。 但 是 在 不 
同 的 参数 和 方案 下 , 相 变 曲线 与 CEP 的 位 置 均 不 尽 
相同 ,在 一 级 相 变 区 域 ,不 同 的 方案 给 出 的 共存 区 大 
小 也 不 尽 一 致 。 但 是 ,这 些 计算 仍然 给 出 了 一 些 确 
定 的 信息 ,首先 , 随 着 化 学 势 增加 , 相 变温 度 减 小 ,这 
证 实 了 前 述 的 由 热力 学 的 一 般 讨 论 给 出 的 结论 ;其 
次 ,所 有 的 计算 均一 致 地 给 出 了 CEP 的 存在 , 进 一 
步 给 出 了 大 化 学 势 区 域 的 共存 区 ,相应 的 结果 展示 
在 图 4 中 。 


0.14 
0.12 
0.10 
% 0.08 
z 0.06 EA 
0.04 | Chiral Symmetry Broken, 


— ~ CEP ~ (0.12, 0.13) GeV 


Chiral Symmetry 
Restored 


0.02 


i \ 
+ Co-eXistence 
1 


(b) 


0 0.1 02 03 04 05 06 
1 GeV 


4 QC 胶 子 模型 下 ,DS 方 程 方法 给 出 的 QCD HIF EAS A RA 
Fig.4 QCD phase diagram under the QC interaction model using the QCD DS equation approach®”! 


目前 ,结合 fRG 方 法 的 DS 方 程 方法 的 最 新 计算 
给 出 了 QCD 的 正确 跑 动 行为 。 首 先 ,在 真空 下 夸克 
质量 函数 随 动量 的 跑 动 行为 与 格 点 QCD 模拟 给 出 
的 结果 一 致 ;其 次 ,在 小 化 学 势 区 域 , 相 变 温度 与 相 
变 曲线 的 曲率 也 有 格 点 QCD 一 致 ,相应 的 结果 已 列 
于 表 1 中 。 因 此 ,可 以 认为 ,目前 的 DS 方程 方法 下 
的 计算 给 出 了 关于 QCD 相 变 的 可 靠 的 结果 。 如 图 5 
BRAN ,相关 计算 再 次 确认 了 CEP 的 存在 ,最 新 的 计 
算 给 出 CEP 状 态 约 位 于 (T,ys)~ (110, 600) MeV. 
2.2” 退 禁闭 相 变 相 图 

再 者 ,DSE 方程 方 法 还 给 出 了 退 禁 闭 相 变相 图 
的 计算 。 在 靠近 相 变 温度 时 ,Schwinger 函数 本 身 表 
PL HH TE ze ME, (ALE Schwinger 函数 的 高 阶 导数 仍然 出 
L 了 负 值 ,最 终 如 图 6 所 示 ,Schwinger 函数 判 据 给 
出 退 禁 闭 相 变 与 手 征 相 变 同 时 发 生 。 男 外 ,通过 
Polyakov 圈 对 胶 子 传播 子 的 计算 也 给 出 了 相同 的 结 
论 。 尽 管 胶 子 部 分 在 Polyakov 圈 中 占据 主要 部 分 ， 
但 是 相 变 的 主要 特征 仍 来 源 于 夸克 质量 函数 的 


kuxi 


S 


变化 。 


et al. | 
0 200 400 600 800 1 000 


5 QCD 的 DS 方程 方法 最 新 给 出 的 QCD 的 手 征 相 变 的 相 


可 
KI 


Fig.5 QCD phase diagram using the DS equation approach in 
the state of the art truncation scheme” 


3 QCD 相 变 理论 和 实验 现象 研究 的 展望 


我 们 知道 ,ALICE/CERN 正在 进行 的 实验 计划 
和 相对 论 重 离子 对 撞 器 (Relativistic Heavy Ion 
Collider, RHIC)/#i 鲁 克 文 国 家 实验 室 
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isoh @ (b) 
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6 QCDIBA AE" (a) Schwinger 函数 判 据 ,(b) Polyakov 势 判 据 
Fig.6 QCD deconfinement phase diagram under the criterion of Schwinger function (a) and Polyakov loop potential (b)""" 

(Brookhaven National Laboratory , BNL) 的 能 量 扫 描 是 重子 数 等 守恒 荷 的 涨 落 及 其 比率 (已 有 初步 结 
以 及 俄 重 离子 同步 加 速 器 (Nuclotron-based Ion 果 52590)、 热 力学 量 和 输 运 系数 ,以 及 夸克 、 胶 子 和 
Collider facility, NICA) #1] HADES/CBM/FAIR 及 中 介子 的 谱 函 数 的 信息 等 。DSE 方 程 方法 在 这 方面 已 
国 强 流 重 离子 加 速 器 (High Intensity Heavy-ion 经 迈 出 的 有 希望 的 一 步 , 并 将 在 未 来 发 展 扩大 ,对 避 
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寻 QCD 相 变 的 CEP 信 号 视 为 其 主要 任务 。 而 如 何 。 致谢 感谢 秦 思 学 教授 、 付 伟 杰 教授 等 的 讨论 和 
JA QCD 的 理论 方法 对 QCD 相 图 的 结构 和 强 相互 作 。 建议 . 
用 物质 的 相关 物理 做 出 可 靠 的 定性 和 定量 预言 是 对 eae, 、 本 
前 理论 研究 的 巨大 挑战 。 其 中 ,关于 QCD 的 DS ”作者 贡献 说 明 两 位 作者 高 飞 、 刘 玉 金 共同 负责 编 
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